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ABSTRACT 

We present a three-component mixing model for the evolution of O 
abundance relative to Fe, taking into account the contributions of the first very 
massive stars (with masses of ^ 100 Mq) formed from Big Bang debris. We show 
that the observations of O and Fe abundances in metal-poor stars in the Galaxy 
by Israelian et al. (1998) and Boesgaard et al. (1999) can be well represented 
both qualitatively and quantitatively by this model. We use the representation 
of the number ratios (0/Fe) versus l/(Fe/H). In this representation, if there is 
only a single source with a fixed production ratio of O to Fe beyond a certain 
point, the subsequent evolution of (0/Fe) is along a straight line segment. 
Under the assumption of an initial Fe ([Fe/H] ~ —3) and O inventory due to 
the prompt production by the first very massive stars, the data of Israelian et 
al. (1998) and Boesgaard et al. (1999) at -3 ^ [Fe/H] ^ -1 are interpreted to 
result from the addition of O and Fe only from type II supernovae (SNII) to 
the prompt inventory. At [Fe/H] ;^ — 1, SNII still contribute O while both SNII 
and type la supernovae contribute Fe. During this later stage, (0/Fe) sharply 
drops off to an asymptotic value of ~ 0.8 (O/Fe)©. The value of (0/Fe) for the 
prompt inventory at [Fe/H] ~ —3 is found to be (0/Fe) ~ 20 (O/Fe)0. This 
result suggests that protogalaxies with low "metaUicities" should exhibit high 
values of (0/Fe). The C/0 ratio produced by the first very massive stars is 
expected to be <C 1 so that all the C should be tied up as CO and that C dust 
and hydrocarbon compounds should be quite rare at epochs corresponding to 
[Fe/H] ^ -3. 

Subject headings: Galaxy: abundances — Galaxy: evolution — supernovae: 
general 



(0/Fe)-Februaxy 1, 2008 



3 



1. Introduction 

In a recent theoretical study of Fe and r-process abundances in the early Galaxy 
by Wasserburg & Qian (2000a), it was proposed that the first stars formed after the 
Big Bang were very massive 100 M©) and promptly enriched the interstellar medium 
(ISM) to [Fe/H] ~ —3, at which metaUicity formation of regular stars (with masses of 
~ 1-60 Mq) took over. Subsequent Fe enrichment was provided initially by a subset of 
Type II supernovae (SNII). The interpretations of Wasserburg & Qian (2000a) were based 
on observations of metal-poor stars in the Galaxy by Gratton & Sneden (1994), McWilliam 
et al. (1995), McWiUiam (1998), and Sneden et al. (1996, 1998). The critical metaUicity 
[Fe/H] ~ —3 deduced by Wasserburg & Qian (2000a) for transition from formation of 
very massive stars to regular stars is in remarkable coincidence with the lower bound 
of [Fe/H] ~ —2.7 observed in damped Lya systems over a wide range in redshift (e.g., 
Prochaska & Wolfe 2000). Furthermore, Wasserburg & Qian (2000b) showed that the 
[Fe/H] of damped Lya systems can be explained in a straightforward manner by addition 
of Fe from SNII to the prompt inventory provided by the first very massive stars and that 
the evolution of [Fe/H] with redshift reflects the chronology of turn-on of protogalaxies. 

The first very massive stars would certainly produce elements other than Fe. One of the 
issues to be resolved is the general problem of the nucleosynthetic signatures of these stars. 
This problem was discussed in the early work by Ezer & Cameron (1971). More recently, 
studies on condensation of Big Bang debris (e.g., Bromm, Coppi, & Larson 1999) indicated 
that protostellar aggregates of ~ 100 M© may be formed. Evolution of zero-metallicity stars 
of ^ 100 M© was discussed by Fryer, Woosley, & Heger (2000). Calculations by Heger, 
Woosley, & Waters (2000) suggest that C, N, O, Mg, Si, and Fe are produced in large 
amounts by these stars. The question is how the theoretical models of nucleosynthesis in 
the first very massive stars formed from Big Bang debris are related to the observations. 
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Another issue concerns the Fe enrichment by Type la supernovae (SNIa). It was assumed 
by Wasserburg & Qian (2000a,b) that SNIa started to contribute Fe at [Fe/H] ~ —1 and 
were responsible for ~ 2/3 of the solar Fe inventory. In a self-consistent Fe enrichment 
model, the history of SNIa Fe contribution needs to be examined more carefully based on 
observations. 

In this paper we present a phenomenological model for the evolution of O abundance 
relative to Fe. Our approach differs from previous studies (e.g., Timmes, Woosley, & 
Weaver 1995) in the following respects. Previous studies rely on many inputs such as the 
nucleosynthetic yields predicted by various models of SNII and SNIa and some description of 
star formation. In the approach here we attempt to identify a chemical evolution model that 
can account for the observed evolution of O abundance relative to Fe in terms of several key 
parameters. We hope that this approach will give a simple yet complete interpretation of 
the data. The usual representation of the O and Fe data is to display [0/Fe] against [Fe/H]. 
In the approach here we use (0/Fe) versus l/(Fe/H), with the round brackets indicating 
the number ratios. If there is only a single source with a fixed production ratio of O to Fe 
beyond a certain point, the subsequent evolution of (0/Fe) as a function of l/(Fe/H) is 
along a straight line segment. As will be shown, this representation also facilitates a simple 
geometric means to describe the evolution of O and Fe abundances in the general case. As 
the first very massive stars are thought to be a source of O and Fe (Heger et al. 2000), 
the value of (0/Fe) at [Fe/H] ^ —3 should serve as a nucleosynthetic signature of these 
stars, thus addressing the first issue raised above. Subsequent to the prompt enrichment by 
the very massive stars, SNll are essentially the only source of O, while Fe enrichment was 
provided by SNII initially and augmented by the contribution from SNIa only at a later 
time. Thus the evolution of O abundance relative to Fe at [Fe/H] > —3 provides a means 
of addressing questions concerning the onset and nature of SNIa Fe contribution. 
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Based on the observations of metal-poor stars in the Galaxy by Israelian, Garcia 
Lopez, & Rebolo (1998) and Boesgaard et al. (1999), the value of (0/Fe) resulting from 
the prompt production by the first very massive stars is ~ 20 (O/Fe)0. Comparison of our 
model with the same observations gives a history of SNIa Fe contribution that is consistent 
with the results of previous studies (e.g., Timmes et al. 1995) and with the assumptions 
of Wasserburg & Qian (2000a,b). The average production of O relative to Fe for SNII is 
found to be (0/Fe)ii ~ 3 (O/Fe)©, which is ~ 7 times lower than that for the very massive 
stars. Our model is presented in §2 and compared with observations in §3. Discussion and 
conclusions are given in §4. 



2. Model for Evolution of O Abundance Relative to Fe 

For simplicity, we focus on the evolution of O and Fe abundances in a closed system 
with a fixed mass. This mass can be stored in gas and stars. The initial state P at time 
t = contains gas only. At t > 0, the total mass of gas Mgas(t) and that of stars Mstar(^) 
satisfy 

Mgas(t) + Mstar(t) =Mgas(0). (1) 

We ignore the return of O and Fe to the gas by stars that are not related to SNII or SNIa 
and consider that stars sample the composition of the gas at the time of their birth. Thus 
the mass fraction X^it) of the element E in the gas is governed by 

XEit)M,Ut) + I XEit') U t' = [PE,ii(t') + P^,m]dt\ (2) 

where -Pe,ii(^) and -Pe,i(^) are the total rates for production of the element by SNII and 
SNIa. 

Differentiating equation (Q) with respect to t gives 

^Mgas(t) = PE,Il(t)+PE,l(t), (3) 
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where equation (|l|) has been used. Equation (|^) can be integrated to give 



XE{t) = Xe(0) + 



Mgas(t') 



dt' + 



Mga.(t') 



dt'. 



(4) 



Note that the integrands in equation only depend on the production rates per unit mass 
of gas. As the mass fraction of H in the gas is nearly constant during its evolution, X^it) is 
proportional to the number ratio (E/H). So equation can be rewritten as 



where (E/H)p is the inventory for the initial state P, and (E/H)n and (E/H)i represent the 
integral of the production rate per unit mass of gas over time for SNII and SNIa. 

We consider that the first very massive stars (with masses of ^ lOOM©) formed from 
Big Bang debris provided an inventory of Fe ([Fe/H] ~ —3) and O. Frequently we use 
"prompt inventory" or "prompt production" when referring to this inventory of Fe and 
associated elements. It is considered that only very massive stars could form until this level 
of Fe and associated elements was reached. The onset of formation of regular stars (with 
masses of ~ 1-60 M0) corresponds to [Fe/H] ~ —3. This defines the initial state P of the 
system. Subsequently, SNII that are related to the formation of regular massive stars (with 
masses of ~ 10-60 Mq) began to contribute O and Fe. At a still later time (see §3), SNIa 
also began to contribute additional Fe. It is considered (e.g., Timmes et al. 1995) that 
SNIa essentially produce no O. So we assume a production ratio (0/Fe)i = for SNIa. 
The O and Fe yields may vary significantly between individual SNII (see §4). However, the 
O or Fe contribution from all SNII may be lumped together so long as we are concerned 
with timescales of, e.g., ^ 10^ yr (see §4). Furthermore, we will assume a constant average 
production ratio (0/Fe)ii for SNII. 

The value of (0/Fe) for the initial state P resulting from the prompt production by 
the first very massive stars is taken to be a constant (0/Fe)p. Thus (0/Fe)p represents 



(E/H) = (E/H)p + (E/H)ij + (E/H) 



I' 



(5) 
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the production ratio for the prompt source, (0/Fe)n for SNII, and (0/Fe)i = for SNIa. 

The whole evolution of O and Fe abundances can then be described by a three-component 

mixing model. In consideration of the time sequence indicated above, there is a domain 

where the model reduces to only two components, the prompt source and SNII (no SNIa). 

In this domain, the number ratio of O to Fe in the gas is 

(n /Pp^ = (Q/H)p + (Q/H)„ ^ (0/Fe)p(Fe/H)p + (0/Fe)„(Fe/H)n 
^""/'"^ (Fe/H)p + (Fe/H)„ (Fe/H) ' 

where (Fe/H) = (Fe/H)p + (Fe/H)jj. Equation (||) can be rewritten as 

(0/Fe) = (0/Fe)ij + [(0/Fe)p - {0/Fe\,]^^^. (7) 

[te/ti) 

Thus the ratio (0/Fe) is a linear function of (Fe/H)p/(Fe/II) with a slope of 
[(0/Fe)p — (0/Fe)jj]. The evolution of (0/Fe) in this regime would be along the 
line segment PII in Figure 1. Point P is at (0/Fe) = (0/Fe)p and (Fe/H)p/(Fe/H) = 1, and 
represents the initial state P. Point II is at (0/Fe) = (0/Fe)ii and (Fe/H)p/(Fe/H) = 0, and 
represents the result after an infinite number of SNII (eq. [0). If SNII cease to occur after 
contributing finite amounts of Fe and O, then the end point for this part of the evolution 
would lie somewhere on the line segment PII. 

If the SNII contributions were to cease and SNIa were to begin to contribute Fe at the 
state corresponding to point A on the line segment PII in Figure 1, then the evolution of 
(0/Fe) would again reduce to a two-component mixing model with one member being state 
A and the other being SNIa. In this simplified case (no concurrent production by SNII), 
the evolution of (0/Fe) beyond point A would be along the line segment AI in Figure 1 
corresponding to 

(0/Fe) = (0/Fe) 



(Fe/H)^ 
(Fe/H) ' 



where (Fe/H) = (Fe/H)^ + (Fe/H)j, and the subscript "A" stands for quantities at point 
A. Point I in Figure 1 is at (0/Fe) = (0/Fe)i = and (Fe/H)p/(Fe/H) = 0, and represents 
the result after an infinite number of SNIa (eq. ^j). 
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As SNII still occur after SNIa start to contribute Fe at point A, the actual evolution 
of (0/Fe) follows a curve between the line segments All and AI in Figure 1. This is a 
three-component mixing regime involving state A, SNII, and SNIa. An important parameter 
for this regime is the rate per unit mass of gas for Fe production by SNIa relative to SNII: 

_ djFe/Rydt 
^ d{Fe/R\,/dt' ^ ' 

For a constant 7, we have 

(0/Fe)„ 



,0/Fe) ^ (2Z^ ^ 

7 + 1 



(0/Fe)^ - 

7 



(Fe/H)^ 



(10) 



(Fe/H) ' 

corresponding to the line segment AC in Figure la. Point C is at (0/Fe) = (0/Fe)n/(7 + 1) 
and (Fe/H)p/(Fe/H) = 0, and represents the result after an infinite number of SNII and 
SNIa (eq. [l^l)- It can be seen that the line segments All and AI correspond to the limits 
7 = and 00. 

In general, 7 is a function of time t. We can discuss the evolution of (0/Fe) as a 
function of (Fe/H)p/(Fe/II) for the case of a time-dependent 7 by calculating d {0/Fe)/ dt 
and d [(Fe/H)p/(Fe/H)]/(it. From the rates for addition to the total O and Fe inventory in 
the gas, 

d (0/H) _ d (0/H)ji _ (o/Fe)ji^^i^^, (11) 



and 



we obtain 



dt dt dt 



d (Fe/H) ^(Fe/H)„ , djFe/R), ^ , ,, d{Fe/}i\, 



_ d{0/Fe) (0/Fe) - {(0/Fe)„/[7(t) + 1]} 

d [(Fe/H)p/(Fe/H)] (Fe/H)p/(Fe/H) " ^ > 



Figure lb gives a simple geometric interpretation of equation ([T3|) . At an arbitrary point 
B beyond point A, the tangent 5* of the curve for (0/Fe) is the same as the slope of the 
line segment BC. Point C is at (0/Fe) = (0/Fe)ii/(7B + 1) and (Fe/H)p/(Fe/H) = 0, and 
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represents the result after an infinite number of SNII and SNIa for a constant 7 equal to 
the value 7b at point B (eq. [0). 

Equation ([1^) also gives 



According to equation (0), if •yit) increases with time, i.e., (i7/(i [(Fe/H)p/(Fe/H)] < 0, 
then (i S/d [(Fe/H)p/(Fe/H)] < and the corresponding curve for (0/Fe) is concave 
downward. On the other hand, if 7(t) decreases with time, i.e., d'-f/d [(Fe/H)p/(Fe/H)] > 0, 
then dS/d [(Fe/H)p/(Fe/H)] > and the corresponding curve for (0/Fe) is concave upward. 
Figure Ic illustrates the trend of (0/Fe) for some hypothetical time evolution of 7(t) 
specified hj d'j /dt = 'y. In the simple case where 7(t) monotonically increases with time, 
the curve for (0/Fe) drops more and more steeply in approaching (Fe/H)p/(Fe/H) = 
(Fig. lb). In any case, the value of (0/Fe) corresponding to (Fe/H)p/(Fe/H) = is 
(0/Fe)n/(7oo + 1), with 700 being the asymptotic value of 7(t) after an infinite amount of 
time. Note that in the above treatment, the evolution of (0/Fe) is directly related to the 
Fe inventory (Fe/H). The time t only serves as an implicit parameter that has the same 
sense of increase as (Fe/H). 



We now proceed to compare the phenomenological model presented in §2 with 
the observational data on the evolution of O abundance relative to Fe. We take 
[Fe/H]p = log (Fe/H)p — log (Fe/H)^ = —3 to correspond to the initial state P resulting 
from the prompt production by the first very massive stars. Here the subscript "©" 
indicates quantities pertaining to the sun. The data from the observations of metal-poor 
stars (mostly unevolved halo stars) in the Galaxy by Israelian et al. (1998) and Boesgaard et 




(14) 



d [(Fe/H)p/(Fe/H)] (Fe/H)p/(Fe/H) [7(t) + 1]^ d [(Fe/H)p/(Fe/H)] " 



3. Comparison with Observations 



(0/Fe)-February 1, 2008 10 

al. (1999) are shown as (O/Fe)/(O/Fe)0 = IQio/Fe] versus (Fe/H)p/(Fe/H) = ioP^/hIp-P^/h] 
in Figure 2a. The data at (Fc/H)p/(Fe/H) < 0.066 are shown without error bars in the 
inset. Note that the data of Boesgaard et al. (1999) have been modified to match the scale 
of stellar parameters used by Israelian et al. (1998) (see Israelian, Garcia Lopez, & Rebolo 
2000). Figure 2a shows that the evolution of (0/Fe) follows a hnear trend between the 
initial state P and an intermediate state at (Fe/H)p/(Fe/H) ~ 0.02. The evolution beyond 
the intermediate state follows a much steeper curve. These evolution trends are expected 
from the model, with the intermediate state corresponding to state A at the onset of SNIa 
Fe contribution. 

For a quantitative comparison of the model trajectory with the data on the evolution 
of O abundance relative to Fe, we assume a time parameterization where the rate per unit 
mass of gas for Fe production by SNII is constant over Galactic history. This rate can be 
written as 

d(Fe/H)„_a(Fe/H)^ ^^^^ 



d t tsSF 

where a is the fraction of the solar Fe inventory contributed by SNII and issF is the time of 
solar system formation (with t — Q corresponding to the initial state P). The occurrence of 
SNIa is considered to depend on the accretion of matter onto a white dwarf from its binary 
companion. The formation of white dwarfs requires long-term stellar evolution and thus 
SNIa cannot occur early in Galactic history. We assume that SNIa Fe production begins to 
occur at state A and then increases to a constant value relative to SNII Fe production over 
a timescale f (it will be shown that the results are not sensitive to the choice of f). This is 
represented by the following expression of the rate per unit mass of gas for Fe production 
by SNIa relative to SNII: 

t~tA 



1 — exp ^- 



r 



^(i-U), (16) 



where is the time corresponding to state A, and 9{t — tjC) is a step function with its 
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value being unity at t > and zero otherwise. Under the above assumptions, the O and 
Fe inventory in the gas is 



(0/H) = (0/Fe)p(Fe/H)p + a(0/Fe)ii(Fe/H)Q(t/tssF) 



(17) 



and 



(Fe/H) 



(Fe/H)p + «(Fe/H)0(t/tssF) 



+ a7oo(Fe/H)p, 



exp 



t-tA 



e{t-tA). (18) 



After choosing a set of parameters (0/Fe)p, (Fe/H)p, (Fe/H)A, a, and t/^ssf, we 
can obtain (0/Fe) as a function of (Fe/H)p/(Fe/H) from equations (P^) and (|TB|). The 
quantities (0/Fe)jj, ^a, and 7oo in these equations are given by 

(Fe/H)p" 



(0/Fe)ji 



1 

a 



tA 



_ (Fe/H)^ - (Fe/H)p 



and 



a7o 



tssF a(Fe/H)0 

1 - - [(Fe/H)p/(Fe/H)^] 



(19) 
(20) 

(21) 



1 - (^aAssf) - (r/tssF){l - exp[-(tssF - tA)/f]}' 
Note again that the time t (measured in units of tssr) is introduced as an imphcit parameter 

to be ehminated from the final results. The form of this parameterization does not 

significantly affect these results. 

As a numerical example, we take (0/Fe)p = 21 (O/Fe)0, (Fe/H)p = 10"^ (Fe/H)©, 
(Fe/H)A = 0.05 (Fe/H)0 [corresponding to (Fe/H)p/(Fe/H)A = 0.02], a = 1/3 for the 
fraction of solar Fe inventory contributed by SNII, and r/tssF = 0.1. The corresponding 
trajectory of (0/Fe) as a function of (Fe/H)p/(Fe/H) is shown as the solid curve in Figures 
2a, 2b, and 2c. These figures cover a series of nested intervals in (Fe/H)p/(Fe/H) to 
expose the complete range of results. For reference, the same results (solid curve) are 
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also shown in the conventional representation of [0/Fe] versus [Fe/H] in Figure 2d. The 
data of Edvardsson et al. (1993) for Galactic disk stars are included in Figures 2c and 2d 
to show the trend near and beyond the solar point (circle). It is evident that the model 
provides both a good qualitative and quantitative description of all the data. There is some 
scatter of the data relative to the solid curve from the model over almost the entire range 
of (Fe/H)p/(Fe/H), especially in the range 0.03 ^ (Fe/H)p/(Fe/H) ^ 0.1. However, we do 
not consider this to be critical. 

To test the sensitivity of the results to the choice of parameters, we have varied 
T^/tssF, (Fe/H)A, or a one at a time while keeping all the other parameters the same 
as in the above example. The long-dashed curve in Figure 2c is obtained by taking 
'T'/tssF = 1 and the short-dashed curve by taking (Fe/H)A = 0.1 (Fe/H)© [corresponding to 
(Fe/H) p/ (Fe/H) A = 0.01]. These two variations seem to have little effect on the description 
of the data. The only significant change is that the value of (O/Fe)/(O/Fe)0 in the infinite 
future decreases from 0.8 for the solid curve to 0.36 for the long-dashed curve due to the 
tenfold increase in f/tssF- Changing a from 1/3 to 1/4 only slightly improves the agreement 
with the data in the range 0.03 ^ (Fe/H)p/(Fe/H) ^ 0.1 (not shown). It would be difficult 
to obtain a fit to the data for a significantly greater than 1/3. 

The value of (0/Fe) for the initial state P appears to be (0/Fe)p ~ 20 (O/Fe)© from 
the existing data. This clearly shows that the contributors to the prompt Fe inventory have 
also made major contribution to O and very possibly to C, Mg, and Si. The value (0/Fe)p 
is considered here to be the result from the production by the first very massive stars 
(with masses of ^ 100 M©) and not by SNII (see §4). The results from complete models of 
nucleosynthesis in these stars should provide a guide as to what may be expected for the 
other elements. The average production of O relative to Fe for SNII is (0/Fe)n ~ 3 (O/Fe)©, 
which is ~ 7 times lower than (0/Fe)p for the very massive stars. The high value of 
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(0/Fe)p required by the data implies that even after the major addition of Fe and O 
from SNII, the total O inventory at the onset of SNIa Fe contribution may still have a 
significant component from the prompt source. For example, at (Fe/H)p/(Fe/H) ~ 0.02 
([Fe/H] ~ —1.3), ~ 10% of the O inventory but only ~ 2% of the Fe inventory are from the 
prompt source. 

We are further concerned with the possibility that prior to the formation of very 
massive stars with masses of ^ 100 Mq, there may have been a separate generation of 
supermassive stars with masses of ^ lO'^-lO'^ Mq. Such objects may produce peculiar values 
of (0/Fe) in the data around (Fe/H)p/(Fe/H) ~ 1. We have sought to find an indication of 
this, but there does not appear to be any such indication in the data. 

4. Discussion and Conclusions 

We have presented a phenomenological model for the evolution of O abundance relative 
to Fe. In this model, the first very massive stars formed from Big Bang debris provided a 
prompt inventory of Fe ([Fe/H]p ~ —3) and O, thus defining an initial state of the ISM. 
Subsequent O enrichment was provided by SNII that are related to regular massive stars, 
while subsequent Fe enrichment was provided by SNII followed by additional contribution 
from SNIa. Prior to the onset of SNIa Fe contribution, the evolution of (0/Fe) as a function 
of l/(Fe/H) is along a straight line segment (Fig. 1). Subsequent evolution is governed by 
the rate per unit mass of gas for Fe production by SNIa relative to SNII (Fig. 1). Based on 
the observations of Israelian et al. (1998) and Boesgaard et al. (1999), the value (0/Fe)p 
resulting from the prompt production by the first very massive stars is ~ 20 (O/Fe)©. 
Comparison of our model with the same observations indicate that SNIa began to contribute 
Fe at [Fe/H] ~ —1 and were responsible for ~ 2/3 of the solar Fe inventory, consistent 
with the results of previous studies (e.g., Timmes et al. 1995) and with the assumptions of 
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Wasserburg & Qian (2000a,b). The average production of O relative to Fe for SNII is found 
to be (0/Fe)ii ~ 3 (O/Fe)©, which is ~ 7 times lower than that for the first very massive 
stars. 

The model developed here neither depends on nor gives numerical values of the rates 
per unit mass of gas for O and Fe production by SNII, d (0/II)jj/(it and d(Fe/}i)^j/dt. 
Likewise, the model only involves 7(t), the rate per unit mass of gas for Fe production by 
SNIa relative to SNII. Taking the Galaxy for example, we can obtain the numerical values 
of (0/H)jj/(it and (i (Fe/H)jj/(it by introducing the time for solar system formation 
tggp ^ 10^° yr (eqs. [|n|] and [T^). The absolute yields of O and Fe from each SNII can 
also be obtained by further introducing the SNII rate per unit mass of gas and the total 
mass of gas at tssp- The progenitor stars of SNII have very short lifetimes so that an 
equilibrium between the birth and death of these stars is established on a correspondingly 
short timescale. If the rate for formation of these stars is proportional to the mass of 
gas, then the SNII rate per unit mass of gas can be considered as constant over Galactic 
history. This is the justification for our earlier assumption of a constant rate per unit mass 
of gas for Fe production by SNII (eq. ||15|). At the present time, the Galactic SNII rate is 
~ (30 yr)~^, corresponding to a total gas mass of ~ 10^° Mq. Assuming that all SNII have 
the same O yield, we find that each SNII must produce ~ 0.3 Mq of O in order to enrich 
the gas with a solar O mass fraction 9.6 x 10"^) at tgsr ~ 10"'^'^ yr (eq. 

In the model presented here, we have assigned the O and Fe abundances at [Fe/H] ~ —3 
to the prompt production by the first very massive stars formed from Big Bang Debris and 
have rejected the possibility that these could be produced by SNII. Indeed, although the 
value of (0/Fe)p found here is ~ 7 times larger than the average relative production of O 
to Fe that we associate with SNII, it is still possible to attribute the large (0/Fe) value 
at [Fe/H] ~ —3 to contributions from a subset of SNII that produce O but little Fe. The 
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key argument that we have used in rejecting this possibihty comes from the results on the 
production of r-process elements and Fe by SNII that were obtained from meteoritic data 
and observations of metal-poor stars. We give a short summary of these results below as 
they will be used to place bounds on the O and Fe enrichment to be expected from the 
early occurrence of SNII in comparison with the prompt inventory (0/H)p and (Fe/H)p. 

Meteoritic data on the inventory of radioactive ^^^I and ^^^Hf in the early solar system 
require at least two distinct kinds of SNII as sources for the r-process (Wasserburg, Busso, 
& Gallino 1996; Qian, Vogel, & Wasserburg 1998; Qian & Wasserburg 2000). These are the 
high-frequency SNIIH events responsible for heavy r-process elements with mass numbers 
A > 130 (e.g., Ba and Eu; "H" for "high-frequency" and "heavy r-process elements") 
and the low-frequency SNIIL events dominantly producing light r-process elements with 
A < 130 (e.g., Ag; "L" for "low-frequency" and "light r-process elements"). The timescales 
required for replenishment of the appropriate radioactive nuclei in an average ISM by the 
SNIIH and SNIIL events are ~ 10^ yr and ~ 10^ yr, respectively. These timescales can be 
explained by considering the mixing of the cjccta from an individual SNII with the ISM. 
The amount of ISM to mix the SNII ejecta is the mass swept up by a supernova remnant 
(SNR). The expansion of an SNR can be roughly described by an energy-conserving phase 
followed by a momentum-conserving one (e.g., Thornton et al. 1998). This expansion 
results in a total swept-up mass of ~ 2£^sN/(i^tri^f) ~ 3 x 10^ Mq, where £^sn ^ 10^^ erg is 
the supernova explosion energy, vtr ~ 200 km s~^ is the velocity at the transition between 
the two phases of expansion, and Vi ~ 15 km is the final velocity typical of the random 
motion in the ISM. For Galactic rates of ~ (30 yr)^^ for H events and ~ (300 yr)~^ for 
L events in a total gas mass of ~ 10^° Mq, we obtain that the rates of SNIIH and SNIIL 
events in an ISM of ~ 3 x 10^ M© are ~ (10^ Y^)'^ and ~ (10^ yr)~^, respectively, as 
required by replenishment of the appropriate radioactive nuclei in this ISM. 
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Observations of metal-poor stars in the Galactic halo (McWiUiam et al. 1995; 
McWilliam 1998; Sneden et al. 1996, 1998) show that there is wide dispersion in the Ba 
abundance loge(Ba) = log(Ba/H) + 12 at [Fe/H] ~ —3 (region b in Fig. 3). This dispersion 
in the Ba abundance (~ 2 dex) at a nearly constant [Fe/H] led Wasserburg & Qian (2000a) 
to conclude that the SNIIH events cannot produce a significant amount of Fe. The Ba 
abundance becomes extremely low for [Fe/H] < —3 and there is no correlation with Fe in 
this region (region a in Fig. 3), indicating that there was some Fe production prior to the 
occurrence of SNII. Observations (Gratton & Sneden 1994; see also McWilliam et al. 1995) 
also show that there is a correlation between the abundances of Ba and Fe at [Fe/H] ^ —2.5 
(region c in Fig. 3). 

Estimates of the Ba abundance resulting from a single SNIIH event (Qian & Wasserburg 
2000; Wasserburg & Qian 2000a) are shown as the band labeled "1 SNIIH" in Figure 3. The 
low Ba abundances in stars with [Fe/H] ~ —4 to —3 (region a in Fig. 3) compared with this 
band and the sharp increase in the Ba abundance at [Fe/H] ~ —3 led Wasserburg & Qian 
(2000a) to conclude that a source other than SNII must exist to produce Fe (along with 
other elements such as C, N, O, Mg, and Si) at [Fe/H] ^ —3. They attributed this source 
to the first very massive stars (with masses of ^ 100 M©) formed from Big Bang debris. 
They further argued that formation of regular stars (with masses of ~ 1-60 Mq) could not 
occur until later when [Fe/H] ~ —3 was reached. In this later regime, regular stars with 
masses of ^ 10 M0 explode as SNII. The more frequent occurrence of SNIIH events with no 
Fe production over a period of ^ 10^ yr gives rise to the sharp increase in the Ba abundance 
at [Fe/H] ~ —3. This is followed by the occurrence of low- frequency Fe-producing SNIIL 
events. Contribution to Fe from SNIa only takes place at [Fe/H] ^ —1. Enrichment in Fe 
of metal-poor stars with —1 ^ [Fe/H] ^ —2.5 must be provided by the SNIIL events. For a 
Galactic rate of ~ (300 yr)~^ corresponding to a total gas mass of ~ 10^° Mq, each SNIIL 
must produce ~ 0.1 Mq of Fe in order to enrich the gas with ~ 1/3 of a solar Fe mass 
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fraction 1.2 x 10"^) at tssF ~ 10^° yr (eq. [Q). Diluting this amount of Fe in an ISM of 
~ 3 X Mq that mixes with the ejecta of a single SNII event results in [Fe/H] ~ —2.5. 
Thus the correlation between the abundances of Ba and Fe at [Fe/H] ^ —2.5 shows that the 
variation of Fe yields between SNIIH and SNIIL events is smoothed out over a timescale of 
10^ yr. 

While the SNIIH events would produce oxygen, they will not be sufficient to provide 
the large O inventory at [Fe/H] ~ —3. Diluting the O yield of ~ 0.3 Mq from an SNIIH 
event into an ISM of ~ 3 x 10"^ Mq would result in (0/H) ~ 10"^ (0/H)q in this mixture. 
We note that two stars with Ba abundances close to our estimate for a single SNIIH 
event and with essentially the same values of [Fe/H] ~ —3 are both observed to have 
(0/H) ^ 10^^ (0/H)q, which is a factor of ~ 10 greater than our estimate for a single 
SNIIH event (HD122563 has loge(Ba) = -1.54 and [Fe/H] = -2.74, and BD -18°5550 
has loge(Ba) = -1.72 and [Fe/H] = -2.91; Westin et al. 2000; McWilliam 1998; Cavallo, 
Pilachowski, & Rebolo 1997). 

We also note that the prompt value (0/Fe)p ~ 20 (O/Fe)© at [Fe/H] ~ —3 assigned 
by us is based on the observations of BD +23°3130 with [Fe/H] = —2.9 (Israelian et al. 
1998), BD +03°740 with [Fe/H] = -3, and BD -13°3442 with [Fe/H] = -3.13 (Boesgaard 
et al. 1999; Israelian et al. 2000) but there are no relevant Ba data. This value of (0/Fe)p 
corresponds to (0/H)p ~ 2 x 10^^ (O/H)©, which could be explained by the O contribution 
from ~ 20 SNIIH events at ~ 0.3 Mq of O per event. If this were the case, then the relevant 
stars should have high log e(Ba) values. Future measurements of Ba abundances in these 
stars are needed to identify possible SNIIH contribution to the observed O abundances and 
to determine the prompt O inventory more accurately. However, from the above arguments 
based on the observations of HD 122563 and BD — 18°5550, we conclude that the origin of 
the O inventory at [Fe/H] ^ —3 must be assigned to the prompt source. 
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Based on the assignment of (0/Fe)p ~ 20 (O/Fe)©, the average production ratio 
of O to Fe for the first very massive stars is ~ 560 by number [(0/Fc)q ^ 28]. As the 
amount of O production is much less than the mass of the star, this means that the typical 
Fe yield from a zero-metallicity star of mass M must be <^ 0.6 M© x (M/IOOM©). In 
considering the prompt enrichment of the ISM by the very massive stars, it is necessary 
to take several effects into account. The universe at the epoch for formation of these stars 
was probably dominated by gas. We do not know how many very massive stars contributed 
to the prompt inventory. Explosions of these stars are likely to be much more energetic 
than supernova explosions. So the mass of ISM for diluting the ejecta from the explosion 
of a very massive star is expected to be ^ 3 x 10"^ Mq. With an O yield of ~ 0.1 M for a 
very massive star of mass M and with the number A^ms such stars contributing to the 
prompt O inventory of (0/H)p ~ 2 x 10~^ (O/H)0, this dilution mass can be estimated as 
~ 5 X 10^ M0(M/1OO Mq) X iVvMS- The expectation that this mass should be > 3 x 10^ Mq 
leads us to conclude that the prompt inventory must have been provided by AVms ^ 1 very 
massive stars. That is, many very massive stars are required to provide (0/H)p. Clearly, 
further studies are needed to understand the formation, evolution, nucleosynthesis, and 
explosion of these stars. 

To determine the relationship between (0/Fe) and (Fe/H) in the Galaxy, Israelian 
et al. (1998) and Boesgaard et al. (1999) used the OH lines of metal-poor stars with a 
wide range of (Fe/H). Their results represent the most complete and self-consistent data 
set available. There is some disagreement between different workers on the evolution of 
O abundance relative to Fe, in part due to the choice of spectral lines used to determine 
the O abundance (e.g., Barbuy 1988). While the general model developed here must stand 
on its own, the numerical value of (0/Fe)p for the prompt inventory is dependent on the 
observational results of Israehan et al. (1998) and Boesgaard et al. (1999). 
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Up to this point, we have assumed that the ISM consists of gas only and that the 
abundances in an individual star are representative of the average gas at the birth of the 
star. In fact, the Galactic ISM consists of both gas and dust, with dust making up ~ 1% 
of the total mass at present. Certain elements such as Fe in the ISM are greatly depleted 
in the gas phase and are considered to be condensed into dust. So long as stars sample 
the composition of the bulk ISM including both gas and dust, we can still interpret stellar 
observations in terms of Galactic chemical evolution as done here because dust is converted 
into gas in stars. In contrast, when applying the model developed here to damped Lya 
systems, it must be further assumed that Fe and O are not substantially fractionated 
relative to each other in the protogalactic ISM of such systems. If for example, Fe were 
dominantly in dust and O in gas in the ISM of damped Lya systems (as is the case in 
many galaxies of about the present epoch) , then a straightforward application of the model 
to such systems would give completely misleading interpretation of the data, which only 
sample elements in the gas phase. 

The most accurate measurement of O abundance in a damped Lya system so far 
was reported by Molaro et al. (2000). They gave [0/H] = —1.85 ± 0.1 corresponding to 
[0/Fe] = 0.19 ± 0.14 for a damped Lya system with [Fe/H] = -2.04 ± 0.1 at a redshift 
z = 3.39. These data do not fit the trend exhibited by the Galactic data shown in Figure 
2d. The data on other elements (e.g., Cr and Zn) in the above damped Lya system also 
suggest that there is no effect of dust depletion (Molaro et al. 2000). We note that the 
observed value of [0/H] in this system is close to the prompt O inventory assigned here 
based on the Galactic data while the value of [Fe/H] is ~ 10 times higher than the prompt 
Fe inventory. We cannot offer a simple explanation for these O and Fe abundances. 

Substantial efforts have been made to establish the dust-to-gas ratio for damped Lya 
systems at redshifts ^ 1. Estimates have been made for the possible dust depletion of 
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condensible elements (e.g., Si, Cr, and Fe) relative to more volatile ones (e.g., Zn) that 
are considered to remain in the gas phase. Pettini et al. (1997) have argued that there 
is some significant dust depletion of condensible elements (by a factor of ~ 2) in the 
damped Lya systems at very early epochs. However, there is strong disagreement between 
various workers (e.g., Lu et al. 1996; Prochaska & Wolfe 2000) regarding the degree of dust 
depletion in these systems. Wasserburg & Qian (2000b) assumed that dust depletion is not 
a dominant effect on the chemical evolution of damped Lya systems at ^ ;^ 1. There is clear 
evidence for large amounts of dust in a variety of galaxies at high z. We cannot reconcile 
the assumption of low dust-to-gas ratios for damped Lya systems with these observations. 
The cause of the transition to the high dust content for galaxies of about the present epoch 
is not understood unless the later dust contribution of Asymptotic Giant Branch (AGB) 
stars plays a major role. 

It is expected that the first very massive stars formed from Big Bang debris produce 
C/0 ratios of <S 1 (Heger et al. 2000). The dominance of O over C should cause all the C to 
be converted into CO, leaving an excess of O in the ISM. This would prevent the formation 
of C grains and also the formation of polyaromatic hydrocarbon (PAH) as well as other 
hydrocarbon molecules. A strict upper limit to the dust content of the ISM at [Fe/H] ~ —3 
can be estimated by taking the total abundance of the major condensible elements Mg, Si, 
and Fe to be ~ 10 times the Fe abundance (see below). As no C dust should be formed, 
dust would constitute at most ~ 10~^ of the total mass of the ISM at [Fe/H] ~ —3. For a 
typical dust size d, the obscuration due to this presence of dust is at most at the level of 
~ 10~^(0.1 iJ,m/d)[N(H. I)/10^^ cm~^] in a system of neutral H column density iV(H I). 

By attributing the O inventory at [Fe/H] ~ —3 to the prompt production by the first 

very massive stars formed from Big Bang debris, the model presented here also predicts that 
the Mg and Si abundances in stars with (0/Fe) ~ (0/Fe)p and [Fe/H] ~ —3 should to a large 
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extent reflect the prompt inventory of Mg and Si. Based on the observations of McWiUiam 
et al. (1995), we obtain (Mg/H)p ~ 3 x IQ-^ (Mg/H)^ and (Si/H)p ~ IQ-^ (Si/H)0, 
corresponding to (Mg/Fe) ~ 3 (Mg/Fe)© and (Si/Fe) ~ 10 (Si/Fe)© for stars with 
[Fe/H] ~ —3. We also consider that the abundances in the prompt inventory represent the 
critical condition for transition from formation of very massive stars to regular ones. Taking 
into account the radiative properties of H, O, Mg, Si, and Fe, we predict that the abundances 
(0/H) ~ 2 X 10-2 (O/H)©, (Mg/H) ~ 3 x lO'^ (Mg/H)©, (Si/H) ~ lO'^ (Si/H)©, and 
(Fe/H) ~ 10""^ (Fe/H)o in the prompt inventory should give sufficient "metals" to permit 
cooling and fragmentation of large condensing gas clouds to form stars with masses of 
1-60 Mq. 

We would hke to thank Garik Israehan for supplying the data of Israehan et al. (1998) 
and Boesgaard et al. (1999) on the same scale of stellar parameters. Comments by the 
referee, F. X. Timmes, were very helpful in improving the presentation. This work was 
supported in part by the Department of Energy under grant DE-FG02-87ER40328 to Y.-Z. 
Q. and by NASA under grant NAG5-4083 to G. J. W., Caltech Division Contribution No. 
8732(1063). 



(0/Fe)-February 1, 2008 22 

REFERENCES 

Barbuy, B. 1988, A&A, 191, 121 

Boesgaard, A. M., King, J. R., Deliyannis, C. R, & Vogt, S. S. 1999, AJ, 117, 492 

Bromm, V., Coppi, R S., & Larson, R. B. 1999, ApJ, 527, L5 

Cavallo, R. M., Pilachowski, C. A., & Rebolo, R. 1997, PASP, 109, 226 

Edvardsson, B., Andersen, J., Gustafsson, B., Lambert, D. L., Nissen, P. E., & Tomkin, J. 
1993, A&A, 275, 101 

Ezer, D., & Cameron, A. G. W. 1971, Ap&SS, 14, 399 

Fryer, C. L., Woosley, S. E., & Heger, A. 2000, ApJ, submitted 

Gratton, R. G., & Sneden, C. 1994, A&A, 287, 927 

Heger, A., Woosley, S. E., & Waters, R. 2000, in The First Stars, ed. A. Weiss, T. Abel, & 
V. Hill (Berlin: Springer), in press 

Israelian, G., Garcia Lopez, R. J., and Rebolo, R. 1998, ApJ, 507, 805 



Israelian, G., Garcia Lopez, R. J., and Rebolo, R. 2000, |astro-ph / 000 1 426 

Lu, L., Sargent, W. L. W., Barlow, T. A., Churchill, C. W., & Vogt, S. S. 1996, ApJS, 107, 
475 

McWilliam A. 1998, AJ, 115, 1640 

McWilliam A., Preston, G. W., Sneden, C, & Searle, L. 1995, AJ, 109, 2757 

Molaro, P., Bonifacio, P., Centurion, M., D'Odorico, S., Vladilo, G., Santin, P., & 
Marcantonio, P. D. 2000, ApJ, submitted 

Pettini, M., King, D. L., Smith, L. J., & Hunstead, R. W. 1997, ApJ, 478, 536 

Prochaska, J. X, & Wolfe, A. M. 2000, ApJ, 533, L5 



(0/Fe)-Februaxy 1, 2008 23 

Qian, Y.-Z., Vogel, P., & Wasserburg, G. J. 1998, ApJ, 494, 285 

Qian, Y.-Z., & Wasserburg, G. J. 2000, Phys. Rep., 333-334, 77 

Sneden, C, Cowan, J. J., Burris, D. L., & Truran, J. W. 1998, ApJ, 496, 235 

Sneden, C., McWilliam, A., Preston, G. W., Cowan, J. J., Burris, D. L., & Armosky, B. J. 
1996, ApJ, 467, 819 

Thornton, K., Gaudlitz, M., Janka, H.-Th., & Steinmetz, M. 1998, ApJ, 500, 95 

Timmes, F. X., Woosley, S. E., & Weaver, T. A. 1995, ApJS, 98, 617 

Wasserburg, G. J., Busso, M., & Gallino, R. 1996, ApJ, 466, L109 

Wasserburg, G. J., & Qian, Y.-Z. 2000a, ApJ, 529, L21 

Wasserburg, G. J., & Qian, Y.-Z. 2000b, ApJ, 538, L99 

Westin, J., Sneden, C, Gustafsson, B., & Cowan, J. J. 2000, ApJ, 530, 783 



This manuscript was prepared with the AAS macros v4.0. 



(0/Fe)-Februaxy 1, 2008 



24 



Fig. 1. — Illustration of the evolution of (0/Fe) as a function of (Fe/H)p/(Fe/H). Here 
(0/Fc)p represents the production ratio for the prompt source, (0/Fe)n for SNII, and 
(0/Fe)i = for SNIa. The initial state P results from the production by the prompt source. 
State A corresponds to the onset of SNIa Fe contribution. The general evolution is described 
by a three-component mixing model involving the initial state P, SNII, and SNIa. Prior to 
reaching state A, the evolution is along the line segment Pll, with point II representing the 
state to be reached after an infinite number of SNII. (a) If SNII contributions were to cease at 
state A, subsequent evolution would be along the line segment PI, with point I representing 
the state to be reached after an infinite number of SNIa. On the other hand, if the rate per 
unit mass of gas for Fe production by SNIa relative to SNII were a constant 7, the evolution 
beyond state A would be along the hue segment AC, with point C representing the state to 
be reached after an infinite number of SNII and SNIa. Note that the hne segments All and 
AI correspond to the limits 7 = and 00. (b) In general, the rate per unit mass of gas for Fe 
production by SNIa relative to SNII is a function of time j{t). So the evolution beyond state 
A follows a curve between the hne segments All and AI. The tangent at an arbitrary point B 
on the curve is the same as the slope of the hne segment BC, with point C representing the 
state to be reached after an infinite number of SNII and SNIa for a constant 7 equal to the 
value 7b at point B. Thus if ^{t) monotonically increases with time, the curve drops more 
and more steeply in approaching (Fe/H)p/(Fe/H) = 0. (c) In general, the curve is concave 
downward if 7(t) increases with time (7 > 0) and is concave upward if 7(t) decreases with 
time (7 < 0). In any case, the value of (0/Fe) to be reached after an infinite number of SNII 
and SNIa is determined by the asymptotic value 700 of "y{t). 

Fig. 2. — Comparison of the model for the evolution of O abundance relative to Fe with 
the data for metal-poor stars in the Galaxy (filled circles: Israehan et al. 1998; squares: 

Boesgaard et al. 1999). Note the change in the trend of the data at (Fe/H)p(Fe/H)~ 0.02 
highlighted in the inset of (a). The data for Galactic disk stars (crosses: Edvardsson et 
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al. 1993) are included in (c) and (d) to indicate the trend near and beyond the solar point 
(circle). The solid curve is obtained from the model by taking (0/Fe)p = 21 (0/Fe)Q, 
(Fe/H)p = 10-3 (Fe/H)0, (Fe/H)A = 0.05 (Fe/H)© [(Fe/H)p/(Fe/H)A = 0.02], a = 1/3, and 
T^/tssF = 0.1. The long-dashed curve in (c) is obtained by taking f/tssF = 1 and the short- 
dashed curve by taking (Fe/H)A = 0.1 (Fe/H)© [(Fe/H)p/(Fe/H)A = 0.01] while keeping all 
the other parameters the same as for the solid curve. 

Fig. 3. — Data (asterisks: Gratton & Sneden 1994; squares: McWilliam et al. 1995; 
McWilliam 1998; triangles: Sneden et al. 1996, 1998) on loge(Ba) versus [Fe/H] for metal- 
poor stars in the Galaxy. The band labeled "1 SNIIH" shows estimates of the loge(Ba) 
value resulting from a single SNIIH event. The wide dispersion in loge(Ba) at [Fe/H] ~ —3 
(region b) indicates that the SNIIH events cannot produce a significant amount of Fe. The 
correlation between the abundances of Ba and Fe at [Fe/H] ^ —2.5 (region c) results from 
the mixture of the SNIIH and the Fe-producing SNIIL events. The low Ba abundances in 
stars with —4 ^ [Fe/H] ^ —3 (region a) compared with the contribution from a single SNIIH 
event and the sharp increase of loge(Ba) at [Fe/H] ~ —3 (region b) suggest that a prompt 
source other than SNII must exist to produce Fe (along with other elements such as C, N, 
O, Mg, and Si) at [Fe/H] ^ -3. 
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